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1 Introduction 

MEG is a useful method to know brain activities, 
and has a remarkable feature that the source location 
is estimated with high accuracy. However MEG is 
insensitive to a radial current dipole, because a 
radial dipole in a sphere makes no magnetic field 
outside the sphere and the shape of human head 
resembles a sphere. To realize more accurate 
localization, there is a requirement to estimate a 
radial current dipole or a dipole with radial 
component in human brain. On the other hand, EEG 
detects the potential at surface of scalp regardless of 
a dipole direction. The accuracy of source 
estimation by using EEG is not good comparing 
with that using MEG, because the human head 
consists of many different electric conductivity 
tissues such as brain, skull and scalp and it is 
difficult to know the electrical conductivity exactly. 
We developed a MEG/EEG hybrid method for 
source localization of a dipole with radial 
component [1]. 

There is a hypothesis that the primary 
somatosensory area 3b and the area 1 respond to the 
median nerve stimulation [2]. The activity in the 
area 1 is considered as a radial current dipole. 

Then, we carried out source location using somato¬ 
sensory evoked field (SEF) and somatosensory 
evoked potential (SEP) by the developed hybrid 
method to verify the hypothesis. The result shows 
that the tangential dipole appears in area 3b at the 
latency of 18-20 ms. The radial component appears 
and increases gradually during the latency of 20- 
23 ms. 

2 Method 

2.1 MEG/EEG hybrid method 

Several MEG/EEG hybrid methods were reported. 
Most of these methods were estimated source 
location by using the weighted data of MEG and 
EEG. However, these methods are only suitable for 
the source estimation of a radial dipole or dipole 
with very small tangential component. Then we 


developed a new MEG/EEG hybrid method by 
focusing the both merits of MEG and EEG. In this 
method, the correlation between MEG data and EEG 
data is adapted to judge whether a dipole has radial 
component. 

The procedure of the developed hybrid method is as 
follows. At first, a sequence of data is divided into 
some sections according to the following items. 

• Case l:Both of MEG and EEG data have 
significant signals in the section. 

• Case 2:MEG data do not have significant signal, 
but only EEG data has significant signal. 

In the period of case 1, the source dipole should 
have tangential component, whether the radial 
component exists or not. We estimate source 
location (x,y,z) and the tangential components 

( m e by using only MEG data. In this step, we 

assume the radial component of the dipole is zero, 
because MEG scarcely include the magnetic field 
due to the radial component. Then, the radial 
component ( m r ) is estimated by using only EEG 
data on the assumption that the location and 
tangential components are the estimated ones by 
MEG. As the results, all parameters of source dipole 
are obtained. 

In the period of case 2, MEG has not significant 
signal, so the source dipole should be a simple radial 
dipole. Then the source location and radial 
component are estimated by using only EEG data. ; 
The number of unknown parameter is 4 in this 
estimation. 

2.2 Measurements of SEF and SEP 

We tried to estimate the response area to median 
nerve stimulation by applying the presented method. 
The SEF and SEP by the stimulation were measured 
with a subject of healthy right-handed and 22 years 
old male in the magnetically shielded room at 
Kohnan-hospital in Sendai, Japan. The electric 
stimulation was a pair of two pulses. The first pulse 
is 0.0-0.2 ms duration and has magnitude of +4 mA, 
and the second pulse is 0.2-0.4 ms duration and has 
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electric current 



Figure 1: Form of stimulation pulse for 
stimulus of the right median nerve 
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Figure 2: Arrangement of electrodes 



Time [ms] 


Time [ms] 


Figure 3: Waveforms of the SEFs using 
Neuromagl22™ induced by stimulus of the right 
median nerve 


Figure 4: Waveforms of the SEPs over the scalp oj 
the contra-lateral hemisphere induced by stimulus 
of the right median nerve 


magnitude of -4 mA. The form of stimulation pulses 
is shown in Fig. 1. This stimulation is delivered at 
the right wrist of the subject with 2.8 Hz inter¬ 
stimulus interval (ISI). The SEFs were gathered 
thorough a low-pass filter of 410 Hz and a high-pass 
filter of 0.03 Hz with 1275 Hz sampling by using 
Neuromagl22™. We recorded SEFs from 50 ms 
before stimulation to 350 ms after stimulation and 
the number of averaged 2000. 

SEPs were recorded by 25 electrodes array 
(5x5 grid of electrodes with 25 mm spacing). We 
placed these electrodes over the scalp of the contra¬ 
lateral hemisphere, so that the center position of 
electrodes array was C3 of the international 10-20 
standard system. The arrangement of electrodes is 
illustrated in Fig. 2. A reference electrode was set at 
the rear of 2.5cm from C4. A low-pass filter of 
410 Hz and a high-pass filter of 0.03 Hz were used 
in the recording. SEP was recorded from 50 ms 
before stimulation to 350 after stimulation and the 
number of averaged 2000. These conditions were 
the identical conditions of SEF. 


SEF and SEP were measured simultaneously. 
Typical to measurements are shown in Fig. 3.and 
Fig. 4. 

2.3 Head model 

In the estimation by MEG, we adopted a single layer 
spherical model. In the estimation by EEG, 3- 
layered meshes were used for head model. These 
meshes were constructed from MR Images taken at 
5 mm intervals. The y-axis connects the nasion and 
the inion, and the x-axis is perpendicular to the y- 
axis at its midpoint and is parallel to the line that 
connects the right and left auditory meati. The z-axis 
is perpendicular to the x- and y-axis through their 
intersection point. 

2.4 Source localization 

We assumed that the source is a single dipole and 
the dipole is in the left hemisphere. On the source 
localization by MEG, we used the combined method 
of the genetic algorithm and the simulated annealing 
method [3]. On the source localization by EEG, we 
used the argumented Lagrange multiplier method. 

























































During the periods between 18 and 23 ms and 
between 30 and 34 ms, the features of SEF and SEP 
were regarded as the easel, while they were 
considered as the case2 during the period between 
24 and 30 ms. 

3 Results 

3.1 Source localization in 18 ms - 20 ms 

Fig. 5(a) shows the result of source localization 
during the period between 18 and 20 ms. We made 
sure by overlay on the subject’s MRI that the dipole 
is on the posterior wall of the central sulcus. The 
dipole is almost tangential dipole and the direction 
of the moment is anteriorly. The change of dipole 
position and the moment vector are little in this 
period. 

3.2 Source localization in 20 ms - 23 ms 

Fig. 5(b) shows the result of source localization 
during the period between 20 and 23 ms. The dipole 
moved gradually to the lateral from the position 
during 20-23 ms. And the radial component 
increased gradually with time, and peaked at 
22.8 ms. This direction change of dipole is listed in 
Table 1. 

3.3 Source localization in 24 ms - 30 ms 

In this period, we could not localize the dipole 
position, because the amplitude of SEP was small. 
Therefore it is supposed that the brain activity is 
very weak in this period. 

3.4 Source localization in 30 ms - 34 ms 

Fig. 5(c) shows the result of source localization 
during 30-34 ms. The dipole is significantly anterior 
(about 4mm) and medially (about 4 mm) from the 
position of dipole during 18-20 ms. The dipole is 
nearly tangential dipole and the direction is opposite 
to that during 18-20 ms. The position and the 
moment vector changed little. 



Figure 5: Source location and orientation of the 
dipole at the latency between: a) 18-20 ms; b) 20- 
23 ms [I- 20.4 ms II- 21.2 ms III- 22.0 ms IV- 22.8 ms] 
and c) 30-34 ms 


T able 1: Position and moment of the dipole at the latency between 20 and 23 ms 




Position [cm] 

Moment [nAm] 


Fatency[ms] 

X 

y 

z 

Radial 

Tangential 

e 

Tangential 

I 

20.4 

-4.95 

1.77 

8.73 

4.83 

-2.67 

-36.79 

II 

21.2 

-5.18 

1.78 

8.68 

14.76 

-2.76 

-27.76 

III 

22.0 

-5.50 

1.72 

8.53 

19.95 

-2.12 

-16.93 

IV 

22.8 

-5.76 

1.85 

8.53 

24.54 

-1.49 

-10.45 














































































4 Discussion 

4.1 Sources in 18-20 ms and 30-34 ms 

A current dipole analysis indicated that the dipole 
during 18-20 ms after stimulus onset was on the 
posterior wall of the central sulcus and oriented in 
anterior. The dipole during 30-34 ms after stimulus 
onset was located medial and superior to the 
position during 18-20 ms and oriented posteriorly. 
These results coincide with the previous MEG 
studies [4,5]. 

4.2 Source in 20-23 ms 

The dipole moved to the lateral and the radial 
component increased gradually from 20 ms to 23 ms. 
From the result, it is speculated that the activated 
area by the median nerve stimulation moves from 
the surface of area 3b to area 1 during 20-23 ms as 
illustrated in Fig. 6(a). However, we assumed the 
source as a single dipole in this work. That is, we 
could not distinguish two dipoles appeared 

simultaneously. Therefore, we regard the source as 
single dipole whose strength, orientation and 

position coincide with the vector addition of two 
current dipoles. Then we could suppose another 
situation as shown in Fig. 6(b). The cortical neurons 
in the area 3b were excited at the latency of 20 ms. 
Subsequently, the cortical neurons in the area 1 were 
excited. The former excitation decreased and the 
later excitation increased with passage of time. 
Though we could not decide which the situation is 
real neurons activity from this study, the presented 
study revealed that radial component of source 
dipole exists between the latency of 20 ms and 
30 ms 
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Figure 6: Illustrations of our speculation of the dipole activety. 

(a) Model l:The dipole moves from area 3b to area 1 with the passage of time. 

(b) Model 2: The cortical neurons in the area 3b were excited. Subsequently, the cortical 
neurons in the area 1 were excited. The former excitation decreased and the later excitation 
increased with the passage of time. The single dipole model approximately leads the position 
of centroid of two dipoles and vector sum of each moment. 
















